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Abstract: 

Many neutron stars (NSs) and runaway stars apparently come from the same regions on the sky. This 
suggests that they share the same birth places, namely associations and clusters of young massive stars. 
To identify NS birth places, we attempt to find NS-runaway pairs that could be former companions 
that were disrupted in a supernova (SN). The remains of recent (<few Myr) nearby (< 150 pc) SNe 
should still be identifiable by observing the emission of rare radioisotopes such as 26 Al and 60 Fe that 
can also be used as additional indicators to confirm a possible SN event. We investigated the origin 
of the isolated NS RX J1605. 3+3249 and found that it was probably born w 100 pc far from Earth 
0.45 Myr ago in the extended Corona- Australis or Octans associations, or in Sco OB4 w 1 kpc 3.5 Myr 
ago. A SN in Octans is supported by the identification of one to two possible former companions - the 
runaway stars HIP 68228 and HIP 89394, as well as the appearance of a feature in the 7 ray emission 
from 26 Al decay at the predicted SN place. Both, the progenitor masses estimated by comparison with 
theoretical 26 Al yields as well as derived from the life time of the progenitor star, are found to be 
« 11M Q . 
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1 Introduction 

There are many young associations and clusters of 
massive stars in the solar vicinity that are potential 
birth places of neutron stars (NSs), hence supernova 
(SN) hosts. The NSs born in those SNe were ejected 
from their parent association or cluster shortly after 
form ation due to a kick in an asymmetric SN explosion 
(e.g. |Burrows fc Hayes|[l996} |Janka fc Mueller||1996 



Janka et a l. 2005; Wang et al. 2006; Kissli nger et al.| 
2009 1. This scenario of NS kicks and ejection from its 



parent association or cluster is supported by the obser- 
vation of large NS pr oper motions that indi c ate high 
space velocities (e.g. Lyne fc Lorimer||1994 Lorimer 



et~aT1[T997| [Hansen fc Phinney||1997| |Cordes fc Cher- 



noff|1998||Arzoumanian et al.|2002||Hobbs et aLf2 005 ) 



For a small numb er of NSs parent associations have 
been suggested (e.g. Hoogcrwerf et al 



2001 



Tetzlaff 



et~aL1[2009l [20T0| |2011a| |Bobylev| |2008| |Bobylev fc 
Bajkova 20091. Due to large uncertainties in the NS 



distances and the unknown radial velocities, the re- 
sults are often not unique ( Tetzlaff et al.|20"T0 2011al. 
Therefore, further indicators are needed to decide on a 
particular birthplace. Such indicators may be the iden- 
tification of a possible fo rmer compan ion that is now a 
so-called runaway star ( Blaauw|1961 that should also 
show signs of the former binary evolution such as high 
helium abundance and high rotational velocity due to 
mass and momentum transfer from the primary as it 
filled its Roche lobe). Other indicators are sources of 
radioactive isotopes. Such isotopes are 26 Al and 60 Fe 



with half-lives of 0.7 2 Myr (e.g. |Rightmire et al.||1958 
Thomas et al.|1984[ ) and 2.62 M^T pugel et al.|2009 | 

respectively, i.e. much longer visible than a SN rem- 



T 



2005 Pons 



tie compar- 



nant (« 10 4 yr). A cooling NS is visible for ss IMyr 
(see e.g. co oling curves in Gusakov et al 
et al.| 2009 1, i.e. a similar time span, 
ison between maps of 7 ray emission, probable ori- 
gins of runaway stars and NSs as well as the SN rate 
shows that there are regions on the sky wh ere more 
SNe/ NSs are present than average (Fig. ll. More- 
over, such regions contain young OB associations, e.g. 
Cygnus, Vela or Orion. Hence, it is plausible to search 
for NS and runaway origins within young associations 
and clusters. Providing small regions on the sky with 
enhanced SN/ NS number is also an important input 
for gravitational wave searches. 

By tracing back in time young NSs, young runaway 
stars and young associations/clusters, we search for 
close encounters between the NS and a runaway star 
and/ or the NS and an association. To account for the 
errors on the observables as well as the for NS unknown 
radial velocity, we perform Monte Carlo simulations. 

If we find that a NS and a runaway star could have 
been at the same place and time inside an association/ 
cluster, it is well possible that the NS was born at that 
place at that time in a SN. From the flight time of the 
NS (its kinematic age), the association age and as- 
suming contemporaneous star formation, we can also 
obtain the mass of the SN progenitor from its life time 
and evolutionary models. Converting the measured 7 
ray flux due to 26 Al decay into a mass of 26 Al that was 
ejected in the SN can then test theoretical nucleosyn- 
thesis yields. 

Here, after describing our method, we present re- 

1 http : / / www .atnf.csi ro.au/research / pulsar / psrcat / 
( [Manchester et al.|2005[ l 
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Figure 1: Top panel: Past flight paths (lMyr 
into the past) of young NSs (red, adopting zero 
radial velocity, kinematic data from ATNF p ulsar 
database 1 ) and m assive runaway stars (blue, 



Tct- 



zlaff et aT]|2011b ). Symbols indicate their present 
position in Galactic coordinates. They are over- 
plotted on t he 26 Al COMPTEL map from |Diehl 
et al. (2010) (contours). Most NSs and runaway 



stars seem to originate from regions on the sky 
with enhanced 7 ray emission (e.g. Cygnus, Vela, 
Orion). Bottom panel: SN rate within 6 00 pc from 
the Sun ( |Schmidt|20lT| see also |Hohle et al.|2010[ ), 
expected from current O and B type stars. 



suits on RX J1605. 3+3249, a member of t he so-called 
'Magnificent Seven" (e.g. |Haberl| |2007| |Mereghetti 



2011 



Kaplan et al. 20111. RXJ1605. 3+3249 is a 



promising candidate for a nearby SN since its distance 
to the Sun is probably < 400 pc ( |Motch et aT1|2007| ). 
Furthermore, for RX J1605. 3+3249 the only way to 
give an estimate on its age (apart from comparison 
with cooling models) is to derive the age kinematically 
since no spin-down wa s measured so far and a possi- 
ble 6.88 s spin period ( jHaberl 20071 is currently not 
confirmed. Precise ages of NSs are also important to 
probe cooling models and constrain their composition, 
hence their equation-of-state. 



2 Procedure 

We followed the same app roach as already descr ibed 
in 



Tetzlaff et al 



(20101 and 



also Hoogerwerf et al 



2001 



Tetzlaff et al. (2011al (see 



Therefore, we limit our 



We apply Monte Carlo simulations varying the pa- 
rameters (parallax, proper motion and radial velocity) 
within their error intervals, to calculate the trajecto- 
ries of the NSs and any association or cluster centre of 
a sample of 296 young associations and clusters within 



3kpc from the Sun into the past [see Tetzlaff et al 



(2011a I for a descripti on of the association sample, for 
associations listed in Mel'Nik & Dambisl (2009), wc 



took their parameters; here we add the Pleiades Bl 



moving group ( Asiain et al. 1999 1 since it was pro- 



posed to have hoste d recent nearby SNe ( jBerghofer fc 



Breitschwerdt 



20021], At every time step (10" up to 



5 x 10^ yr in steps of 10 4 yr) the separation A between 
the association/cluster centre and the NS is calculated. 
We then find the minimum separation A m i„ (r) and 
the associated time r in the past. For the radial ve- 
locity of the NS, we assume a reasonable probability 
distribution derived from the pulsar space velocity dis- 
tribution by |Hobbs et al.| ( |2005[ ). 

The distribution of separations Amin is supposed 
to obey a distribution of absolute differences of two 
3D Gaussians with (mi,si) and (rri2,S2) being their 
expectation values and standard deviations (i f the po- 



sitional errors we re Gaussian distributed), see Hooger- 
werf et al. 



(2001), 



W 3D (A) 



A 



1 (A- 
2 



1 (A_ 
2 



(1) 

where A denotes the 3D separation between two ob- 
jects (here, NS and association centre or NS and run- 
away star; A = A m in), m = |mi — TO2 and s 2 — 
si + sh yr = 3.1459.... 

For investigating encounters with runaway stars, we 
calculate A m i„ between the NS and the runaway star. 
If the two stars once were at the same place, i.e. 
m — > 0, |cquation l| becomes 



W 3 D,m^0 (A) 



2A 2 
2^s 3 



exp 



A^i 

'2s 2 



(2) 



Note that even in this case, there is zero probability of 
finding a value of A m in = in the Monte Carlo sim- 
ulation. However, near-zero values should be found 



description here to the most important facts. 



after a sufficient number of runs (see appendix A I. 
The peak of the distribution is shifted towards larger 
A m ; n depending upon s. As the actual (observed) case 
is different from this simple model (no 3D Gaussian 
distributed positions, due to e.g. the Gaussian dis- 
tributed parallax that goes into the position recipro- 
cally, complicated radial velocity distribution, etc.), we 
will fit the theoretical formulae only to the first part of 
the A min distribution (up to the peak plus a few more 
bins) such that the slope and peak can be explained 
(in the following, we use the term "adapt" instead of 
"fit" because [equations 1| and [2] are only used to ex- 
plain the first part of the A m in distribution, hence it 
is not a real fit). The parameter m then gives the po- 
sitional difference between the two objects. Note that 
the uncertainties on the separation are dominated by 
the kinematic uncertainties of the NS that are typically 
of the order of a few hundred km s _ (because of the as- 
sumed radial velocity distribution). As a consequence, 
the distribution of separations Amin shows typically 
a large tail for larger separations. However, the first 
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part of the A m ; n distributions (slope and peak) can 
still be explained well with |equation 1| since the kine- 
matic dispersions for only those runs are much smaller, 
typically a few tens of kms -1 for the NS, i.e. a few 
tens of pc after 1 Myr. 

To associate the NS/runaway star encounter posi- 
tion with an association/cluster, the trajectory of the 
association/cluster is calculate d simultaneously. Run- 
away star data are taken from |Tetzlaff et aL (2011b I 
(and references therein, mainly [van Leeuwen|2007[ ) for 
2547 runaway stars (1705 with full 3D kinematics). 

In general, we first perform 10 4 Monte Carlo runs 
for each NS/association pair and 10 3 Monte Carlo runs 
for each NS/runaway star pair (for the latter less runs 
are still sufficient due to the smaller errors on the run- 
away star kinematics compared to the dispersion of the 
association velocities) to find those associations and 
runaway stars that potentially crossed the past path 
of the NS, i.e. those for which the smallest A m in value 
found in the calculations is less than three times the as- 
sociation radius (for NS/association pairs) or less than 
10 pc (for NS/ runaway star pairs), respectively (see 
appendix A} . Those associations and runaway stars 
that fulfilled those conditions are then selected for a 
more detailed investigation (one to three million Monte 
Carlo runs). The outcome of these simulations is then 
discussed in detail. Regarding associations, we search 
for those for which the NS could have been within the 
association boundaries in the past while for runaway 
stars, we are looking for those runaway stars for which 
the NS and the runaway star might have been at the 
same place in the past, hence the distribution (slope) 
should obey |equation 2| After three million runs the 
smallest A m j„ value found is exp ected to be smaller 
than one parsec (see appendix A[ ). If this criterion is 
satisfied, we adapt [equations l| and[2|to the first part of 
the "observed" A m i„ distribution to explain its slope 
and peak. If we find that a distribution with m — 



( equation 2 1 can satisfactorily explain the slope and 
peak of the A m i n histogram, we consider the runaway 
star as former companion candidate. 

This procedure was already successfully applied 
by |Hoogerwerf et al.| |2001|), |Bobylev| fl2008|>; |Bobylev| 
fc Bajkova| ( |2009| ) and us QTetzlaff et al.||2009| |20Toj 
2011a[) and was also applied to several (artificial) test 
cases ( |Tetzlaff|200"9 l. 



3 RXJ1605. 3+3249 

We adopt the following parameters 
for RX J1605. 3+3249 (right ascension a, declination 
S, distance d, proper motion fi% — pi a cos 8, us): 



16 h 05 m 18!5, S = + 32°49'17"4 
( Kaplan et al.|2003D, 
350±50pc flPosselt et al.|2007|f 



/'a 

Us 



43.7 ± 1.7 mas/yr jZane et ai.|2006[ ) 
148.7 ± 2.6 mas/yr ( |Zane et al.|2006[ ). 



First, we perform 10 4 Monte Carlo runs to find 
close encounters between RX J1605. 3+3249 and any 
association/cluster in the past five million years 3 ! 
We select those associations/clusters for which the 
smallest separation A min found was less than three 
times the association/cluster radius R aS s, 18 in total. 
For those 18 associations/clusters, we carry out an- 
other one million Monte Carlo runs. For 10 of them 
we find close encounters consistent with the associ- 
ation/cluster boundaries. They are listed in |Table l[ 
We adapt [equation l] to the first bins of each A m i n dis- 
tribution to obtain the distance of the SN to the asso- 
ciation centre. Comparing the radii of each association 
with this putative separation of the SN from the asso- 
ciation centre (m ± s), three associations/clusters are 
found to be potential birth places of RX J1605. 3+3249, 
i.e. the position of the SN is consistent with the asso- 
ciation boundaries within its standard deviation: the 
extended Corona- Australis association (Ext. R CrA), 
Scorpius OB4 and the Octans association. In |Tablc~2"1 
we give the position of the SN and the properties 
RX J1605. 3+3249 would currently have if it was born 
in the respective association. In the last Column we 
give an estimate of the mass of the progenitor star 
derived from the life-time of the progenit or star [as- 
sociation age (Ext. R CrA: 10 — 15 Myr [N euhauscr 
|et al.||2000"l Oct ans: sa 20 Myr, |Torres et al.||2008T 



Sco OB4: 7Myr |Kharchenko et a l. 2005) mi nus N S 
age (t)1 and evolutionary mode ls from Tinsley ( 1980 ) 
Maeder fc Meynetj ( |1989[ ) and |Kodama| ( |1997| . We 



regard these three associations as likely birthplaces of 
RXJ1605. 3+3249. The derived uncertainties in the 



parameters ( Table 2 1 are dominated by the uncertain- 



ties in the NS radial velocity and distance. They are 
included in the formal errors (68%). For more distant 
associations (>few hundred pc) such as Sco OB4, the 
results are also influenced by the uncertainties in the 
association distances. These mostly affect the derived 
position of the SN within the association. Note that 
although the proposed progenitor mass of 42 — 89 Mq 
is rather high (because the lifetime of the progenitor 
star is short assuming it formed in Sco OB4 at the 
same time the association formed) and black holes are 
expected to form above ~ 25 Mq ( Heger et al. 2003 1 



rather than NSs, for masses higher than ~ 50 — 80 Mq 



in binary systems, NSs are expected to form (Belczyn- 
|ski fc Taam|2008[ ). 

After 10 3 Monte Carlo runs 40 runaway stars with 



(3) 



- ^Posselt et al.|2007| used different models for the hydro- 
gen col umn density and derived values of 390 pc and 325 pc. 
IMotch et al.|(|2007l give an upper limit of 410 pc. According 
to |Posselt et al.|2007| the different models give consistent, 
hence very reliable results up to £3 270 pc (their values are 
also in good agreement with the parallactic distances for the 
two NSs RXJ1856. 5-3754 and RX J0720. 4-3125). Hence, 
our adopted distance for RX J1605. 3+3249 of 350 ± 50 pc 
is reasonable. 

3 T he effective temperature of RX J1605. 3+3249 o f kT = 
96 eV jMotch et aTll999| |yan Kerkwijk et al |2004[ | imply 
a cooling age of 10 —10" yr (see e .gl cooling curves by 
|Gusakov et al.|2005| |Pons et al.|2009"] | , hence older ages can 
be excluded. 
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Table 2: Potential parent associations of RXJ1605. 3+3249. 



Assoc. 


(m, s) 




v r 


v-i 




dist 


d 




s 






H 


[Myr] 


[kma- 1 ] [mae/yr] 


[mas/yr] 


[kms- 1 ] 


[PC] 


[PC] 


[°] 


[°] 


[M ] 


Ext. R CrA 
Sco OB4 a 
Octans 


(33.4, 2.3) 
(42.7, 18.4) 
(103.5, 25.7) 


42 + 07 
-0.04 

, 4 + 0.3 

3 - -0.6 

53+ - 09 


577+^ 3 -43.9±1.7 
-21 + Jj> -42.3±1.3 
548 + 3^ 9 -43.6±1.7 


148.4±2.6 
149.8±2.5 
148.6±2.6 


612 + lf 
266+H 


303 + |° 
385 + |° 


ioi± 7 12 

983 + 3 ° 


261.5+l'j 

a «T..+g : ; 


-44.6+2-J 

-34 3+°" 9 
34.3_ g 

-78 ... - 68 


12 - 18 
42 - 89 
10 - 11 


(m, s) — expectation value a 


nd standard dc 


Nation of the encounter 


cparation i 


iforrod from 


equation l| 


putative SN 


distance from a 


ssociation cent 




[m ± a); t - 




Predicted pre 


sent NS parameters: rad 


al velocity 


j r . proper n 




nd fj,g , current distance dist; 


v sp ~ predictc 


d 



space velocity v sp . Predicted SN position: SN d istance d^ (at the time of the SN), J2000.0 coordinates (as seen from Earth at present). Error bars 
denote 68 per cent confidence (cf. Appendix B in |Tetzlaff et al.|2010[ . M prog - progenitor mass 

a Thc results for Sco OB4 were o btained u sing a uniform radial velocity distribution for the NS instead of the one derived from the |Hobbs et al.|(2005[ 
distribution, sec also comment in [Table T] 



Table 1: Properties of potential close encounters 
between RXJ1605. 3+3249 and the centres of 10 
associations/clusters. 



Assoc./ 


Russ 


(m, s) 


T 


cluster 


[pc] 


[pc] 


[Myr] 


Tuc-Hor 


50 


(130.4,65.5) 


63 +UUU 

u -°°-0.11 


P Pic-Cap 


57 


(135.2,59.0) 


60+ 07 

UDU -0.12 


ext. R CrA 


31 


(33.4,2.3) 


42+°-° 7 

U -* Z -0.04 


AB Dor 


43 


(155.0,51.1) 


o =^+°-° 6 

uoo -0.10 


Sco OB4 a 


30 


(42.7, 18.4) 


3 4+°' 3 


Columba 


85 


m > 200 


ss 0.6 


Carina 


58 


(129.2,60.0) 


68+ 11 

u - uo -0.12 


Octans 


111 


(103.5,25.7) 


>53+ ° 9 
U - OJ -0.08 


Argus 


74 


(148.7,52.6) 


uoo -0.09 


Pleiades Bl 


54 


(63.0,4.1) 





(m, s) — expectation value and st andard dev iation of the 
encounter separation inferred from [equation *T| putative SN 
distance from association centre (m ± s)\ r — time of the 
encounter in the past. 

a Using the |Hobbs et al.| (2005j distribution for NS space ve- 
locities in the case or Sco UB4 we found that the absolute 
radial velocity v r of RX J1605. 3+3249 must be very small 
if Sco OB4 was the birth association of RX J1605. 3+3249. 
To achieve better statistics and a clear peak in the A m i n 
distribution to be able to adapt the theoretical curves, we 
repeated the calculations assuming a uniform v r distribu- 
tion. The results for the latter are given here. 



full kinematics from Tetzlaff et al.| (201 lb I show a 
smallest separation A m i„ between RX J1605. 3+3249 
and the runaway star of less than 10 pc. For them, 
another three million runs are performed. We select 
those stars for which separations smaller than one par- 
sec are found after three million runs, 28 in total. To 
further reduce the number of potential former compan- 
ion ca ndidates, we follow the idea of |Chmyreva et al"| 
(20101. They used a reference probability as signif- 



icance indicator that two objects approach to small 
distances taking into account the spatial density dis- 
tribution of the objects (we refer to details on this 



method to Chmyreva et al. 20101. Chmyreva et al 



( 2010 1 investigated pairs of pulsars whereas we inves- 
tigate NS-runaway star pairs. Therefore, beside the 
spatial density distribution of NS, also the density dis- 
tribution of runaway stars is needed. As our region 
of interest lies in the Solar neighbourhood, the pul- 



sar distribution used by Chmyreva et al. (20101 is not 



sufficient for our purposes as it does not account for 
an increased SN rate in the Solar neighbourhood and 
includes all NSs (from young to very old). For that 
reason, we use the spatial distributions of young NSs 
and runaway stars derived by Tetzlaff (20121. They 



away stars are ejected from their parent associations 
and clusters. As parent associations and clusters, an 
observed sample of associations and clusters was taken 
and member stars distributed according to the initial 

~ ( [2005 1 



mass function from Kroupa & Weidner 



SN 

times were then obtained from the age of the associa- 
tions and the progenitor lifetime (for details we refer 
to Tetzlaff 20121. From the present spatial positions 



of simulated NSs and runaway stars, their spatial dis- 
tributions were derived. 

From the Monte Carlo simulations we obtain the 
region {X + AX, Y + AY, Z + AZ) in which the NS 
and runaway star can reach separations A m i n no larger 
than p. If p <§; A (X,Y, Z), the number of objects 
in the volume p 3 obeys a Poisson distribution. Thus, 
the spatial region (X + AX, Y + AY, Z + AZ) is di- 
vide d into N cells (fc = 1 . . . k) with edge lengths p 
(see Chmyreva et al.]|2010 l. Using the density distri- 
bution of runaway stars and NSs derived earlier, the 
number of stars (t)ns, Vrun) in each cell is calculated 
and an average number in a cell of size p is obtained. 
The probability of detecting N objects of population i 
(NSs or runaway stars) in a volume p 3 is given by 



Po (p) = E 



(4) 



The probability of detecting at least one NS and one 
runaway star in a volume with size p 3 is then given by 

P (p) = (l-e-"« s ) (1-e-"— ). (5) 

The probability of occurrence of separations A m in < p 
in the Monte Carlo simulations is P(p). An upper 
limit Pn.upp of Po can be derived from the Poissonian 
error on the calculated number of stars in each cell. If 
Po,u P p < P for small p (smaller than a few pc), this is 
an indicator that the association of RX J1605. 3+3249 
with a particular runaway star might be real. 

Using this as a criterion, none of the associations 
between RXJ1605. 3+3249 and the former 28 com- 
panion candidates is found to be highly significant. 
However, for those cases, where the encounters are 
found within an association, the reference distribu- 



performed a population synthesis where NSs and run- 



tion for NSs and runaway stars (Tetzlaff 20121 pre- 
dicts a relatively high number of stars, as expected. 
Hence, those encounters are expected to be found 
less significant. Therefore, we still regard runaway 
stars for which encounters with RX J1605. 3+3249 are 
found within one of the potential parent associations 
of RX J1605. 3+3249 listed in |Table 2| namely the Ext. 
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R CrA, Sco OB4 or the Octans association, as good 
candidates]^] Here, for the Ext. R CrA and Octans 
associations, we do not only regard the nominal radii 
(maximum extensions) of the associations but their ex- 
tensions in different directions in space taking into ac- 
count the motion of the association (for extensions see 
Fernandez et al. 120081 for Ext. R CrA and lTorres et al.l 



81696, we may exclude it as a companion candidate 
to RXJ1605. 3+3249. 



2008 for Octans). The five remaining former compan- 
ion candidates to RX J1605. 3+3249 are discussed in 
detail in the following paragraphs. For all of them, 
the Octans association is the proposed host associa- 
tion of the SN. This is mainly due to the large size of 
the association^ It also still could be that more dis- 
tant runaway stars that could confirm a SN scenario 
for RXJ1605. 3+3249 in S co OB4 are not present in 
the Tetzlaff et al. (2011b I catalogue for runaway star 



candidates as it is based on Hipparcos data (complete 
to ~ 0.5 kpc). We stress that the non-identification of 
a possible former companion for RX J1605. 3+3249 for 
a SN in Sco OB4 could be caused by the lack of more 
distant runaway stars in the runaway star catalogue or 
by the possibility that the progenitor star was a single 
star. 



HIP 57787 is a metal-poor G5 giant star with a 
metallicity of [F e/H]= -2.10 + 0.16 ( |Anthony-Twarog 
|fc Tw arog 1994). Due to the low metallicity, the run- 
away nature of this star is probably not due to a SN 
event in a close binary system. 



HIP 80448 is a G2V pre-main sequence binary 
with a companion that is itself a single-lined spectro- 
scopic binary ( |Cutispoto et a l. 2003) with an age of 
17.2 ± 2.9 Myr (ITetzlaff et al.||2011bb. Rs age is con- 



sistent with the age of the Octans association (20 Myr, 
Torres et aL|2008| ) which could have hosted the SN. 



HIP 81696 is an 06.5V double star. An encounter 
of the system with RX J1605. 3+3249 could have oc- 
curred inside the Octans association if the present 
distance of the runaway star candidate was 36tj| pc. 
However, distance estimates from Ca-II column den- 
sities indicate a distance of the system of ~ 1 kpc. 
Also, the Hipparcos parallax of 17.95+ 16.68 mas (van 



Leeuwen 20071 that we used in our simulations is 
very uncertain. Due to the large distance of HIP 

4 Note, that there might have been associations or clus- 
ters in the past that are already dispersed. However, since 
RX J1605. 3+3249 is very young, its parent association or 
cluster should still exist as the dispersion time scale of as- 
sociations and clusters is much larger (tens to hundreds of 
Myr). 

5 Although |Torres et al.| | |2008| question whether the Oc- 
tans association is expanding, it could well be that it was 
smaller in the past, i.e. at the time of the SN. However, 
we neglect the expansion of associations here to not miss 
candidate stars. The chance of finding former companion 
candidates is naturally higher for large associations, how- 
ever we regard each association with a runaway star indi- 
vidually. Every companion star is discussed in detail and 
further observation of the runaway star might be necessary 
to confirm or reject a particular SN scenario. 



HIP 82977 is an AO (Malkov et al. 20061 eclips 



ing binary of Algol type and lies below the theoretical 
zero-age main sequence (ZAMS) of the models used by 
|Tetzlaff et al.| ( |201 lb[ ) to determine stellar ages, it was 
shifted towards the ZAMS and treated as ZAMS star. 
Therefore, it's age is very uncertain. Since AO stars 
spend a few hundred Myrs on the main sequence, it 
is unlikely that HIP 82997 is only a few tens of Myrs 
old. Hence, it is probably not the former companion 
of RXJ1605. 3+3249. 

HIP 86228 is a variable star with spectral type 
Fill ( see |Sowell et afl|2007) . |Q'Mara fc Simpson| 
( |1972[ ) give a helium abundance for HIP 86228 of 
Nhb/Nh ~ 0.1 that is comparable to the mean he- 
lium abundance of B type stars (for a review see also 
Lyubimkov|2010" I . However, they used a spectral type 
of B1III, hence a much different effective temperature 
for this star, the given helium abundance might thus 
be wrong. Re-measuring the helium abundance of HIP 
86228 might hence be crucial. 

Considering the individual discussion above, we re- 
gard HIP 80448 and HIP 86228 as possible fo rmer co m- 
panion candidates to RX J1605. 3+3249. In |Fig. 2| we 
show the distributions of minimum separations A min 
and corresponding flight times r (the kinematic age 
of the NS) for those runs where both, the NS and 
the runaway star were within the boundaries of the 
Octans association. We adapted the theoretically ex- 
pected distribution ( [equations 1| and [5| to the first 
part of the histogram. For HIP 86228 the curve pre- 
dicts that the runaway star could have been at the 
same place at the same time in the past as the NS, 
whereas for HIP 80448 a fly-by is suggested. Hence 
HIP 86228 is so far our best former companion candi- 
date for RXJ1605. 3+3249. 

Since there are also 842 runaway star candidates 
listed in the Tetzlaff et al. (2011b I catalogue without 



radial velocities available in the literature, we search 
for possible encounters with RX J1605. 3+3249 among 
these stars assuming v r = ±500kms _1 for the run- 
away star. For 47 of them it is possible to find a 
past position as close as 10 pc to RX J1605. 3+3249 af- 
ter 10 4 Monte Carlo runs. After further three million 
runs for those stars, we found 26 candidates with a 
smallest separation to RX J1605. 3+3249 less than one 
parsec and a peculiar space velocity smaller than at 
least 180 km s -1 (that is « 6a above the maximum 
of the distrib ution of runaway star velocities, Tetzlaff 
et al. 2011b I. For those stars, we evaluate the ref- 
erence probability PcQ (see above). We find no sig- 
nificant association between a runaway star candidate 
without known radial velocity and RX J1605. 3+3249. 



e At this point, we adopt a velocity dist ribution for these 
stars according to |Tetzlaff et aC] ( |2011b| |. Otherwise, the 
probability of occurrence of values A m i n , P, is a priori 
small and cannot be compared to Pq. In the case of NSs, 
the large variety of possible velocities is already accounted 
for. 
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(a) HIP 80448 




10 12 14 16 18 20 



(b) HIP 86228 



i" 0.5 





Figure 2: Distributions of minimum separations 
A min and corresponding flight times t for encoun- 
ters between RXJ1605. 3+3249 and the runaway 
star candidates HIP 80448 (a) and HIP 86228 
(b). The solid curves drawn in the A min his- 
tograms (bottom panels) rep resent the th coreti- 
cally expected distributions (equations 1 or [2]), 
adapted to the first part of each histogram, see 
Isection 21 For HIP 80448 we find that there could 
have been a close fly-by (to = 4.7 pc, s = 1.8 pc) 
at 0.33^q q2 Myr in the past (solid line in the r his- 
togram, left panels), however no direct encounter. 
For HIP 86228 the adapted curve predicts that 
RX J1605. 3+3249 and the runaway star could have 
been at the same place (to = Opc, s — 2.3 pc) at 
0.45^0 06 in the past. Both encounters would 
have occurred inside the Octans association. 



Table 3: Properties of potential encounters inside 
Octans between RXJ1605. 3+3249 and our former 
companion candidates. 



HIP to [pc] s [pc] 



[Myr] 



10137 
83003 
89394 
92218 
95800 



77.0 


5.3 




24.9 
8.9 
2.3 
1.3 
2.0 



u - 4D -0.07 

40+ 06 
n 47+0-03 

n 44+0-04 

U '^-0.06 

n 43+0.05 



1: Hipparcos identifier 
2,3: parameters of equations 1 



or 



Col. 
Cols. 

[2] adapted to the first part of the A r 
distributions 
Col. 4: flight time of the stars (kinematic 
age of RXJ1605. 3+3249) 



give the parameters of the adapted curves (m, s) and 
the corresponding flight times r. For three runaway 
stars we find that the NS and runaway star could have 
been at the same place and time within the Octans 
association. 

HIP 83003 is listed as an O type Cepheid variable 
(Simbad databas^]). However, its high V band mag- 
nitude V — 10.62 mag suggests that the star might 
be much more distant than its parallactic distance 
of 163^37 pc (infers absolute magnitude of Mv = 
4.6 mag; for an O type star, My ss to —5 mag, hence 
its distance is w 1 — 13kpc). 

HIP 98500 is an F2 main sequence star. As it 
lies just below the theoretical model ZAMS, hence was 
treated as ZAMS star by iTetzlaff et all d2011bb. Tak- 



ing into account the uncertainty in distance, we es- 
timate its age to be roughly 600^4qo Myr (with ages 
possible as small as w 50 Myr) from its position in the 
HRD. Hence, it is probably too old to be the possi- 
ble former companion of RX J1605. 3+3249. Further- 



more, HIP 98500 is a metal- weak star (Bidelman & 
McConnell|1978 1, suggesting that, if it is a true young 



runaway star, it was ejected due to dynamical interac- 
tions in a dense cluster rather than a SN in a binary 
system. 

HIP 89394 was classified as FPU star b y |Houk fc| 
|Cowley||i"975| >, and revised by |01sen| ( fl979| > to be Am, 
i.e. it is metal-rich. We suggest further investigation of 
that star and treat it as a possible former companion 
candidate. 



We still regard stars for which possible encounters with 
RX J1605. 3+3249 are found with in one of the poten- 
tial birth associations of the NS ( Table 2 1 (for them, 



it is expected to find encounters less significant, as 
explained earlier), as former companion candidates. 
For them, we evaluate whether |equation 2| could fit 
the first part of the histogram of separations A m i n 
between RX J1605. 3+3249 and the runaway star (for 
those Monte Carlo runs yielding small A m i„ values 
within the potential parent association, s ee also|Fig. 2| 
for examples of such distributions). In Table 3 we 



4 Summary and conclusions 

Many NSs and runaway stars apparently come from 
the same regions on the sky, i.e. regions where mas- 
sive star associations and clusters are present (Fig. 1 1. 



To identify birth places of NSs, we attempt to find 
NS-runaway star pairs that could have been former 
companions that were disrupted during the SN of the 
primary. 



7 http://simbad.u-strasbg.fr/simbad, operated at CDS 



www.pu blish . csiro. au/jo urnals /pasa 
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Table 4: Predicted current parameters of 
RX J1605.3+3249 and SN position and time. 



former comp. cand. 
v r [kms -1 ] 
current NS dist. [pc] 
Mq [mas yr" 1 ] 
fi s [mas yr" 1 ] 
v ap [kms -1 ] 

dist. of SN to the Sun [pc] 
Gal. long. £ [deg] 
Gal. lat. b [deg] 
time in the past [Myr] 



HIP 86228 

303±g 
-43.7 ± 1.6 

148.4 ±2.6 

-186 
-75 
92 4 



626Z 
303^ 



587Z 



, + 15 

342.9+1;! 

-1.3 
-0.8 

0.45 ± 0.06 



-7.3"* 



HIP 89394 

;+209 

-56 
,+54 
-40 

-43.7 ± 1.7 
148.6 ±2.6 

-183 
-84 
,+ 16 
-20 
;+2.7 
'-1.0 
,+2.7 
-1.6 



657Z 
99j 
335.1 
-16.7! 
0.43 ± 0.05 



Rows 1-5: NS parameters, Rows 6-9: predicted 
SN position and time. HIP 89394 would have a 
radial velocity of v r nip = 117j 



7+H km s" 



Here, we presented our results on the origin of the 
isolated NS RX J1605.3±3249. We suggest that it was 
born in the Octans association « 0.45 Myr ago. Such 
a scenario is supported by the identification of two for- 
mer companion candidates that are now runaway stars, 
namely HIP 86228 and HIP 89394. For the latter, no 
radial velocity has yet been published. 

For both cases, i.e. either HIP 86228 or HIP 89394 
being the former companion of RX J1605. 3+3249, the 
predicted current parameters of RX J1605. 3+3249 and 
the predicted SN position (coordinates as seen from 
Earth at presen t) and time for this ca se are summa- 
rized in Table 4 (see Tetzlaff et al. 2010 for their deriva- 
tion). 

As the predicted SNe would have been very recent 
(« 0.45 Myr ago) and nearby (« 100 pc), we expect to 



find 7 ray emission from ^°A1 decay (e.g. Diehl et al. 
20101). In the case of HIP 86228, at the predicted po- 



sition the emission is dominated by emission from the 
galactic inner region. Also in the case of HIP 89394, 
the predicted SN position is not far from the galac- 
tic p lane; however, in the COMPTEL 1.8 MeV map 
(e.g. Diehl et al. 2010 I there is a feature centred at 
(I, b) = (334"5,-16"6). According to SN remnant ex- 
pansion theory (Sedov- Taylor expansion, Snowplough 
expansion, IClarke & Carswelll 120031), a SN remnant 
expanding into the interstellar medium with a volume 
density of n = 1 cm -3 would have a size of « 46 pc af- 
ter 0.43 Myr. At a distance of 100 pc, this corresponds 
to an angular size of the SN remnant of ~ 25°. The 
size of the feature is close to the COMPTEL resolu- 
tion (ldeg 2 ). If we treat it as a point source, we ob- 
tain a total flux of 2.5 ■ 10 -5 photon s cm" 1 s~ 1 . Wit h 



a half-life of 2o Al of 0.72 Myr (e.g. Rightmire et al 



1958 Thomas et al.|1984 ), this yields an ejected mass 
oP^Al in the SN of 2.2 • 1Q- 5 M^. Compared to the- 



oretical Al yields by Woosley fc Weaver ( 1995 1 and 



Limongi & Chieffi (20051, this corresponds to a mass 
of the progenitor star of « 1 1 M© . 

The progenitor mass can also be derived from the 
life-time of the progenitor star which can be estimated 
as the difference between the age of the parent associ- 
ation and the age of the NS. With the Octans associ- 
ation being the parent association of the NS and run- 
away star, we estimate a progenitor mass of 10 — 11 Mq 
for both cases. This is also in excellent agreement with 




mass [M J 



Figure 3: Present mass function of the Octans as- 
sociation (black circles) in comparison with the 
initial mass function £ (IMF, solid line) as given 
by Kroupa & Weidner ( |2005[ ), their equation 2. 
Dashed and dotted lines represent the lcr and 1.5a 
IMF boundaries. Within 1.5a we expect up to one 
star with 10 M Q ; the probability for higher mass 
stars is low. 



the estimates from 26 Al for the case of HIP 89394. A 
progenitor mass of 10 — 11 Mq corresponds to spectral 
type Bl on the main sequence. All high probability 
Octans members in the list by Torres et al. ( 2008 1 are 



F to K type stars, hence very low mass stars. From 
the Octans present mass function it might be possible 
within 1.5cr that there was one 10 Mq star in Octans 



(Fig. 31 



Our first scenario involving HIP 86228 as former 
companion candidate might be supported by the exis- 
tence of the old (few 10 4 to 10 5 yr) nearby (< 2 pc) 
SN remnan t RCW 114 ( [Bedford et~afl [l984| [Klin 
eTallpoTo] at (l,b) « (334", -6") (|Green||2009[~A~ 



Poghosyan, priv. comm.). 

We plan to take higher resolution spectra of HIP 86228 
and HIP 89394 to measure their (He) abundances, ra- 
dial velocities and rotational velocities. 
We note that more distant runaway stars that could 
support a SN in Sco OB4 might have been not found 
because of the lack of more distant runaway stars in 
the runaway star catalogue which will be updated and 
extended to more distant stars in the future. 



It has been suggested by Ellis et al. ( 1996 1 that 



a nearby SN should deposit radioactive isotopes such 
as 60 Fe on Earth. Since our proposed SN would have 
been close to Earth (~ 100 pc), we may expect to find 
a signal of 60 Fe. However, there is no significant signal 
(younger than 1 Myr ) detec ted i n the sediment i nvesti - 
gated by|Knie et aL] (|2004[) and|Fitoussi et al.|([2008|. 



This may has several reasons: as the ejecta material 
reached the Earth, it contaminated other parts of the 
surface than the regions where the sediment investi- 
gated by IKnie et all (120041) and IFitoussi et all (12008 ) 



was found, maybe where the conditions were not suit- 
able to allow such sedimentations. The sediment could 
be destroyed over the time or it takes more time for 
the ejecta be up-taken by the ocean and/ or the host 
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material than currently expected. Another possibil- 
ity is that the amount of 60 Fe that is released in the 
SN is much less than assumed or the SN is asymmet- 
ric. Hence, a non-detection of an 60 Fe Peak does not 
exclude a nearby SN event at the particular time. 

If an identification of the birth place of a NS can 
be made, the kinematic age gives the best estimate 
on the true NS age since the spin-down age gives only 
an up per limit in most cases (see also Tetzlaff et al. 
2010| . Moreover, for RX J1605.3+3249 no spin-down 



has been measured yet and the age can only be de- 
termined kinematically. Our result of Tkin ~ 0.45 Myr 
is in good agreement with cooling models of NSs (see 
|Tetzlaff et aLl[20To| for a high magnetic field consist- 
ing of polar an d toroidal compo nents, cooling curves 
fit even better, Pons et al.|2009 1. 
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A On choosing A min limits 

We construct a case such that the NS (assuming zero 
radial velocity for RX J1605. 3+3249) and a star with 
typical runaway star parameters occupied the same 
place 1 Myr in the past and test how small A m i„ can 
be achieved in the Monte Carlo simulation applying 
typical errors on the kinematic parameters. We in- 
vestigate the NS and a simulated runaway star with 
typical runaway star kinematics as well as one with 
an assumed radial velocity of ± 500kms _1 , i.e. our 
adopted range when v r is not known. |Table 5| gives 
the adopted U, V and W values for the runaway stars 
in each case. 

Already after 10 3 — 10 4 Monte Carlo runs separa- 



tions below 10 pc are reached (Fig. 4 1, below lpc after 
« 10 6 runs. To be conservative we carried out many 
more runs (10 6 — 3 x 10 6 ) to not miss runaway stars 
near the limits. 
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Tabic 5: Adopted heliocentric velocities for simu- 
lated runaway stars for four different cases. 



Case 


U [kins- 1 ] 


V [kms" 1 ] 


W [kms" 1 ] 


A 


-11±5 


-10±5 


-5+5 


B 


-11±435 


-10+100 


-5±225 



A - RX J1605. 3+3249 and a runaway star 
at 425^25 P c with fully known 3D kine- 
matics, B - RX J1605.3+3249 and a run- 
away star at 425^^5 P c with an adopted 
radial velocity of ± 500 km s -1 . 




10° 10 1 10 z 10 3 1 4 1 5 1 1 

number of runs 




10° 10' 10 z 10 3 1 4 1 5 1 1 



Figure 4: Evolution of A min during a Monte Carlo 
simulation for constructed cases where the NS and 
the runaway star were at the same place 1 Myr 
in the past. Dashed lines mark the standard de- 
viation. Top panel: Case A - RXJ1605. 3+3249 
and a runaway star at 4251^5 P c with fully 
known 3D kinematics. Bottom panel: Case 
B - RX J1605. 3+3249 and a runaway star at 
425^25 P c w ^h an adopted radial velocity of ± 
500 km s" 1 . 



